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IN SITU PLY STRENGTHS: AN INITIAL ASSESSMENT 


C. C. Chamis* and T. L. Sullivan** 

NASA Lewis Research Center 
Cleveland. Ohio 44135 

ABSTRACT 

The in situ ply strengths in several composites were calculated using a computa- 
tional procedure developed for this purpose. Laminate fracture data for appropriate 
low modulus and hi^ modulus fiber composites. were used in the laminate analysis 
in conjunction with the method of least squares. The laminate fracture data were 
obtained from tests on Modmor-I graphite/ epoxy, AS-gpraphite/epoxy. boron/^epoxj’ 
and E- glass/ epoxy. The results obtained show that the calculated in situ ply strengths 
can be considerably different from those m asured in unidirectional composites, 
especially the transverse strengths and those in angleplied laminates with transply 
cracks. 


Key words: fiber composites, graphite composites, boron composites, glass 
composites, fracture stresses, uniaxi.il data, in situ ply strength, stress analysis, 
laminate analysis, least squares method, computer program. 
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INTRODUCTION 

It has been suspected for some time now in the fiber composites community 
that uniaidal ply strength in ang^eplied laminates (in situ ply strengths) may be 
consider 2 d>ly different from those measured in unidirectional composites under 
uniaxial loadii^ (uniaxial ply streng&s). A major reason for this suspicion is 
that the predicted fracture stresses of ang^eplied laminates, based on uniaxial 
ply strengths, are generally considerably lower &an measured data. To obtain 
better agreement between theory and experiment, investigators have either 
aibitrarily increased the transverse and shear strengths by some factor (ref. 1), 
modified them in some manner (ref. 2), or neglected them completely (refs. 3 
and 4). A procedure, therefore, is needed which can be used to assess the in 
situ ply strengths and compare these with uniaxial values. To meet this need, 
an investigation was conducted at NASA Lewis Research Center with the objective 
of obtaining an initial assessment of in situ ply strengths. 

The investigation included both experimental and theoretical aspects with 
emphasis on developing a formal procedure to determine the in situ ply strengths 
indirectly. The experimental program consisted of testing quasi- isotropic 
an^eplied laminates from Modmor-l/epoi^, at various load angles. The laminate 
configurations were: (0/±60) (0/ ±45/90) (0/±30/±60/90) and(0/±22.5/ 

S 5 S 

±45/±67.5/90) . Each of these laminates was tested at load angles of 0^, 15^, 

30**, 45°, 60°, 75°. and 90° to the 0°-ply direction. The major portion of this 
work was reported in reference 2. The theoretical part is described in detail 
herein. The analytical effort consisted of using laminate analysis in conjunc- 
tion with the method of least squares to determine the in situ ply strengths using 
fracture stress data from the above quasi- isotropic composites and data available 
In the literature for AS-graphite/epoxy, boron/epoxy and E-^ass/epoxy. 

THEORETICAL APPROACH AND GOVERNING EQUATIONS 

The approach used to calculate the in situ ply strengths is as follows. It is 
assumed that the ply combined- stress strength criterion is given by the following 
complete quadratic function of the in- plane stresses. 
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where the (T| 's denote ply stresses at that load condiUon, the A's denote uniaxial 
strength coefficients which are initially unknown and are to be determined as will 
be described later. The subscripts refer to material axes with 1 taken along the 
fiber direction. It is important to note that equation (1) is not biased. All its 
terms appear with equal weight. Note also that equation (I) permits for coupling 
(interaction) between the normal stresses and <’'£29 and the intralaminar 

shear stress ^^*^ation (1) has been referred to in the literature as the 

strength tensor (ref. 5). 

The A's in equation (1) are determined using the following method of least 
squares matrix equation (ref. 6). 

(<^ll*^l<^£l {A} =((X£l '*'{!} (2) 

9xm mx9 9x1 9xm mxl 

where the <jj^s in equation (2) are the ply stresses at laminate fracture and are 
computed using laminate theory, and where m is greater than 9 and denotes the 
number of identical plies for which distinct ply stresses have been computed. 

By making the following substitutions 

rxf==(<r£)^(0£l (3) 

9x9 9xm mx9 

{y} =l<^£l’^{l) 14) 

9x1 9xm mxl 

equation (2) can be written thus 

(x]{A}={y} (5) 

where the dimensions have been dropped for convenience. 

The A's are determined from equation (5) as follows: 


{A} =(xl’^{y', 


(6) 
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where [ x]~^ may be determined using any of C'le standard matrix inverters such 
as Gauss elimination for example. Once the A's are known, die ply uniaxial 
fracture stresses (in situ ply strengths) are determined by solving equation (1) 
for the various special cases (determine its roots), that is by assuming that the 
ply is subjected to one stress only. For the case where and ^ 

* 0, equation (1) reduces to 




( 8 ) 

( 9 ) 

( 10 ) 


COMPOSITE SYSTEMS EXAMINED 

The composite systems examined were: E-glass/epoxy (E>G/Epon S26 or 
E-G/E (ref. 7)), hi^ modulus graphite/epoxy (Modmor-l/ERL.A 4617 or MOD- 
I/E (ref. 2)), boron/epoxy (boron/ epoxy AVCO 5505/4 or B/E (ref. 8)), and low- 
modulus graphite/ epoxy (AS/3501 or .A$/E (refs. 3 and 9)). The unidirectional 
composite (ply) properties of these composite systems are summarized in table I. 
The reasons for selecting these composite systems are as follows (refer to 
table I): 

1. The E-G/E was selected for its low longitudinal modulus (E^ = 44. l GPa 

(6.4 Msi)), its relatively low transverse tensile strength = 22.7 MPa 

(3. 6 ksi)), and its relatively low intralaminar shear strength * 19-3 MPa 

(1. 5 ksi)). 
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2. The MOD-I/E was selected for Its high (241 GPa (34. 9 Msl)). Us 

relatively low (28 MPa (4.0 kslU. and the hif^ ttiermal expansion coeffi- 

cients difference (A^ = 48x10*^ K”^ (27xi0~® F~^)) which tends to produce transply 
cracks in angleplied laminates when cooling down from cure temperature to room 
temperature. 

3. The B/E was selected for its high (201 GPa (29.2 Msi». Its relatively 

high longitudinal compressive strength = 1592 MPa (231 ksi» and Its reUv 

tively high (55.8 MPa (8.1 ksl)). 

4. The .4S/E was selected for its Intermediate E^^^ (131 GPa (19.1 Msl)>, 
its low transverse modulus (Ej.,o = 8.6 GPa (1.25 Msi», and its relatively high 

S,22t 

5. The extensive amount of laminate fracture stress data available for these 
composite systems (as will be described below),, whieh is required to perform .. 
meaningful least squares analysis, was also a consideration. 

The laminate fracture stress data for the AS^E. B/E and E-G/E composite 
systems were obtained from the literature (refs, cited above) while the data for 
MOD-l/E were generated in-house and reported in reference 2. The various 
laminate configurations and their corresponding loading conditions are summarized 
in table 11. Note that in this table, the distinct ply stress states at fracture for 
each laminate and the total ply stress states for the laminates are given in the last 
t>i-o columns. The number of ply stress states for each laminate is obtained by 
multiplying the number of different plies in each laminate by the number of loading 
conditions for that lamioate. In the case of E-(VE, there were several replicates 
per test and several tests at combined- loading conditions. In table n only the 
totals are shown. The Interested rentier is referred to reference 7 for all the 
details. .As can be seen In table II, a large number of distinct ply stress states is 
avaiLable to evaluate the unknown coefficients (.A's) in equation (1). Since there 
are 9 .A's in equation (1), there are about 3 pi.c strt states per .A for the B/E com- 
posite system ;md about IT for the MOD-L^E. 
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DETERMINATION OF IN SITU PLY STRENGTHS 

The in situ ply strengths for Qie four composite systems (AS/E, B/E, E-G/E, 
and MOD-I/E) were determined using die following procedure. 

1. The various ply stress states in table II were calculated using the linear 

laminate analysis available in the Multilayered Filamentary Composite .Analysis 
(MFCA) Computer Code (ref. lO). The use of linear theory for calculating the 
ply stresses in these laminates is justified because their stress-strain curves 
to fracture show linear behavior, or very nearly so (refs. 2, 3, 7, and 8). The 
inputs required in MFCA are: laminate configuration, laminate fracture stress, 
ply elastic properties ^£l‘> ‘'ll**’ thermal ex- 
pansion coefficients and tablet), and the temperature difference 

between cure and room temperatures (about 167 K and (300® F) for epoxy com- 
posites). 

2. The ply stresses calculated in item (1) above for each composite system 

were used to generate the elements in the [ ] array in equation (2) . For 

example, the first row in [ ] for the .A£/E composite system is the ply stress 

at fracture in the 0®- ply when the laminate was loaded to fracture along the fiber 
o ^ 

(0 -ply) direction, or Note that each row in ( (T^ ) will have 1, 2, or 

3 non-zero elements corresponding, respectively, to ply stresses 
and The other rows in { o-^ ) are generated in a similar manner. The 

dimensions for the array [ Oj j are (60x9) for the A^E composite system (25x9) 
for B/E, (142x9) for E-G/E and (152x9) for MOD-l/E where the first dimension 
((m) eq. (2)) is the total number of ply stress-states corresponding to the last 
column in table n. 

3. The array ( <x^ \ was transformed to generate the elements in the ( x) and 
in the X y } arrays using equations (3) and (4) , 

4. The coefficients .4 were evaluated fi-um equation (6) using Gauss elimin- 
ation. 

5. The desired in situ ply strengths were determined using equations (8), 

(9), and (10) for and ^(22C’ ^il2S’ respectively. 
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The computer program flow chart for ttie above procedure is illustrated 
schematically in figure 1. It is important to keep in mind that the in situ ply 
strengths calculated using this procedure are the best as determined by the 
method of least squares, where ply fracture stresses from uniaxial and angle- 
plied laminates were used with equal wei^t. Equation (1) can be used to study 
the degree of interaction between the normal stresses and the shear stress and 
other features once the A's have been determined. This is discussed in some 
detail in the section General Comments and Recommendations. 

RESULTS .^ND DISCUSSION 

The relative significance of the coefficients Ag and (eq. 1)) which couple 
the normal ply stresses and with tiie intralamihar shear ply stress 

and Ag, which brings in the linear influence of assessed by 

solving equation (6) both with and without these coefficients. The results for 
the AS/E angleplied laminates are summarized in table III. .As can be seen, 
the remaining coefficients are unchanged and, therefore, are independent of 
the presence of .A-, Ag. and .Ag. bimilar results were obtained for the cor- 
responding coefficients of die other composite systems. The significant con- 
clusions from these results are: U) there is no interaction betiveen the in 
situ ply stresses or 0-^22 ‘^£12’ (2) the linear part of does 

not contribute to the in situ ply combined- stress failure criterion; and (3) the 
six coefficients (.Aj^, .A^, A^, A^, .Ag (eq. 1>) are sufficient to describe the in 
situ ply combined-stress failure criterion. In view of this, only these 6 coef- 
ficients were used to determine the in situ ply strengths for all four composite 
system s- 

The predicted in situ ply streng'‘’is for the four composite systems are 
summarized in table IV. The corresponding strengths measured from uniaxial 
tests are also summarized in adjacent columns and the ratios of in situ to 
uniaxial itrengths are given between columns of table IV for comparison pur- 
poses. Note ll»at the in situ ply transverse strengths for the MOD-I/E are not 
given. The reason for this is that the lamination transverse residual stresses 



8 


in tiiese angleplied laminates were of sufficient magnitude to cause transply 
cracks ^ ctescribed in the Composite Systems Examined section. 

Therefore, the plies could not resist additional stress in the transverse direction 
and were not permitted to do so in the laminate analysis during the ply fracture 
stress calculations. 

The following are observed from the results in table IV. 

(1) The predicted in situ ply longitudinal tensile strength Is within 

10 percent of the uniaxial test value for the .AS/E. B/E and E-G/E composite 
systems; however, that for the MOD-I/E system is 22-percent greater. 

2. The predicted in situ ply longitudinal compressive strength is 

within about 10- percent of the unixxial value for the B/E and E-G/E composite 
systems, 12-percent greater for the .AS/E system, and 35-percent smaller for 
the MOD-I/E system. 

3. The predicted in situ ply transverse tensile strength j) is considerably 
greater than the uniaxial test data, about 3 times for .A6/E. 6 times for B/E and 

7 times for E-G/E. As was mentioned earlier, that for the MOD-I/E system 
was taken to be zero because of the transply crack due to high lamination residual 
stresses. 

4. The predicted in situ ply transverse compressive strength is 18-percent 

less than the uniaxial test data for AS/E, 123-percent higher for the B/E system 

and 19-percent higher for the E-G/E sj'stem. That for MOD-I/E was taken to 
be zero for the same reasons as for the transverse tensile strength above. 

5. The predicted in situ ply intralaminar shear strengths aro the 

same as the uniaxial for .AS/E. 1-percent less for B/E, 3-.times greater for E-G/E. 
and 64 percent less for MOD-I/E. 

A graphical comparison of predicted off-axis tensile strength, using uniaxial 
and in situ strengths in the combined- stress failure criterion described in ref- 
erence 11, is shown in figure 2 for all four composite systems. As can be seen 
the difference between the two off-axis strengths is very significant at higher 
load angles. 
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Several important conclusiona maybe made from the above observations. 

1. The in situ ply longitudinal tensile and compressive strengths and 

^IllC ^ taken to be the same as the uniaxial test data for angleplied laminates 
which are free of transply cracks as a result of high residual or Initial stresses. 

2. The in situ ply transverse tensile strength S^22T several times 

higher than the uniaxial data in angleplied laminates free of transply cracks. It 
follows from diis conclusion that calculated laminate strength based on firsts ply- 
failure criteria can be ultra conservative. 

3. The in situ ply transverse compressive strength S^22C either less 

than, greater than or about equal to the uniaxial test data in angleplied laminates 
free of transply cracks. Perhaps uniaxial test methods used to determine this 
strength need re-examination. 

4. The in situ ply shear strengdi be taken to be the same as the 

uniaxial for high- modulus composite angleplied laminates which are free of transply 
cracks. However, for E-G/E it can be several times higher than the unia.xial 
value. 

5. Advanced composite angleplied laminates with trrj'.spiy cracks will tend 

to have a significantly larger in situ ply longitudinal tensile strength than 

measured in uniaxial tensile teats, and considerably lower and Sjj^2s* 

6. Angleplied laminate designs based on first- ply failure and uniaxial ply 
strengths will be conservative if no transply cracks are present; they will also be 
conservative in the presence of transply cracks if the dominant fracture mode is 
controlled by longitudinal tension; however, these laminates will be unconserva- 
tive if transply cracks are present and the dominant fracture mode is controlled 
by either longitudinal compre.<«sion, or intralaminar shear, or both. 

7. The measured values of the uniaxial ply strengths ^{22C 

1‘^S sensitive to the test method used to determine them, and considerable 

care may be required to determine a correct (reasonable) value for these strengths. 
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GENERAL COMMENTS AND RECOMMENDATIONS 

The present investigation led to several significant results and important con- 
clusions concerning in situ ply strengths. From these, some general comments 
can be made which may help shed some light on fiber composite anglepiied laminate 
behavior and at the same time stimulate further discussion and research in this 
area. 

The procedure described herein is general enough and can be used when deter- 
mining in situ ply fatigue life, environmental degradation and not sensitivity. It 
may also be used to assess the adequacy of present uniaxial strength test data, as- 
suming that the in situ ply strengths ar*' more representative of composite behavior 
in structural components. 

It is generally agreed in the composites community that transply cracks have 
no effect on fiber- dominant anglepiied laminates subjected to monotonic or cyclic 
tensile loads. However, the results obtained herein show that transply cracks 
can cause severe reduction in in situ ply longitudinal compression strength and 
intralaminar shear strength. These reductions may be one of the reasons why ad- 
vanced fiber composite anglepiied laminates tend to exhibit poor compression 
fatigue resistance (ref. 12) and severe degradation in residual compress ’-'e strength 
of damaged or flawed anglepiied laminates (ref. 13). The increase in the in situ 
ply longitudinal strength observed herein, in the presence of transply 

cracks, may contribute to the increase in residual tensile strength exhibited by 
advanced anglepiied laminates which have been subjected to tensile fatigue iref. 14). 

Though some aspects of the in situ ply strengths such as interaction between 
normal and shear stresses were investigated herein, several still remain to be 
e.xamined. Some of these are: the goodness of fit of the combined- stress failure 
function, the properties of the solid described by the function iprlnclpal xxes of 
ellipsoid for e.xample), the degree of interaction between normal stresses 
and grraphical comparisons of the function obtained herein with those of 

available failure criteria (similar to fig. for example), sensitivity of the results 
obtained usiair -vailable failure criteria and in situ ply strengths as compared 
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with results obtained using uniaxial strengths, and comparisons of predicted laminate 
fracture stresses with in situ ply strengths versus measured data. 

The procedure described herein can be used to determine the interaction be- 
tween inplane and out-of-plane ply strengths. For this, out-of-plane stresses 
produced from flex, short beam shear, and Interlaminar stresses at free edges 
are needed in conjunction with in plane stresses predicted by laminate throy. The 
in situ ply strengths for other composite systems such as Thomel-300 graphite/ 
epoxy, Kevlar/ epoxy, S-^ass/epoxy end hybrids should be determined using this 
procedure. Finally the possibility of using this procedure to characterize the 
unidirectional composite should be given serious consideration. 

CONCLUSIONS 

The most significant conclusions of an investigation to obtain an initial assess- 
ment of in situ ply strengths are summarized below. 

1. The calculated in situ ply strengths can be considerably di^erent from 
those measured in uniaxial tests, especially those for transverse strength (as 
much as 7 times greater) and intralaminar shear strength (as much as 3 times 
greater). 

2. The calculated in situ ply compressive and intralaminar shear strengths 
can be considerably lower (about 30 and 60 percent, respectively) in angleplied 
laminates with transply cracks; however, the longitudinal tensile strength can be 
higher (about 20 percent) . 

3. A computational procedure consisting of laminate analysis in conjunction 
with the method of least squares can be used to calculate the In situ ply strengths 
and the degree of interaction between these strengths. This procedure is general 
and should be equally applicable for ush in assessing ocher in situ ply strength 
properties such as fatigue, creep, relaxation, notch sensitivity and environmental 
degradation. 

4. The uniaxial ply transverse tensile and compressive strengths and intra- 
laminar shear strength may be sensitive to the test methods used. The current 
practice for measuring these strengths using uniaxial tests is not representative 
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of the In situ state. Perhaps a better alternative may be to adopt calculated in situ 
ply strengths as uniaxial strengths for use in failure load calculations for composite 
components. 

3. The predicted fracture stresses of angleplied laminates, designed using 
uniaxial strengths, will be conservative if they a»;< free of transpiy cracks and 
unconservative if they have transpiy cracks and arc subjected to compressive loads. 

6. The in- plane ply normal strength and intralaminar ply shear strength are 
not coup'ed in the in situ state. 
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TABLE I. - TYPICAL PROPERTIES OF TOE UNIBIRECTKBIAL CCBiPOSTES <PLY) USED IN THIS INVECTIGATION 


Property 

symbol 

Units 
9 <cusQ 

Composite | 

Glasa/epoxy 
E-^asa/epon §26 

High- modulus 
graphits/epoxy 
Modmor l/erla 4617 

(MOD^l/E) 

Boroo/epoxy 
Aveo SWV4 

(B/E) 

Low modulus 
graphli«/«pMy 
W3501 
(AVE) 

Density 

<’1 

gm/cm^ (Ib/ln^l 

0.216 (0.077) 

0.163 (0.056) 

0.205 (0.073) 

0.160 (0.OS7) 

lllMluillS 








®iii 

GPa (Msi) 

44.1 (6.4) 

241 (34.9) 

201 (29.2) 

131 (19.1) 

Transverse 

^i22 


13.8 (2.0) 

7.7 (1.12) 

22.1 (3.2) 

8.60 (1.25) 

Shear 


»» 

4.1 (0.6) 

6.1 (0.89) 

5.4 (0.78) 

3.73 (0.54) 

Poiaaoa*8 ratio 







Ms|or 

‘'112 

— 

0.32 

0.27 

0.17 

0.30 

Thermal exp. coef. 







Longitudinal 

“'ll! 


3.8 (2.1) 

-0.9 (-0.5) 

6.1 (3.4) 

-0.56 (-0.31) 

Transverse 

“|22 

♦# 

36 (20.0) 

46.8 (26) 

30.4 (16.9) 

2.3 (0.13) 

Strength 







Long, tensile 

^niT 

MPa (ksi) 

827 (120) 

563 (81.7) 

1371 (199) 

1630 (236) 

Long, compres. 

®(11C 

»» 

356 (51.7) 

456 (66.2) 

1592 (231) 

1050 (152) 

Transv. tensile 

^22T 


10.3 (1.5) 

28 (4.0) 

55.8 (8.1) 

48 (7) 

Transv. comp. 

®I22C 

ft 

62.7 (9.1) 

200 (29.0) 

128 (17.9) 

255 (36.9) 

Intral. shear 

®I128 

»» 

22.7 0.6) 

51 (8.0) 

62.7 (9.1) 

66 (9.6) 

Ply thickness 


cm <in.l 

0.021 (0.0084) 

0.019 (0.0075) 

(0.005) 

0.014 (0.0055) 


















^WGINAD 
POOR QUAinS 


TABLE n. - SraMARY OF PLY STRESS STATES USED DJ THE ANALYSIS 


Composite system sod 
laminate conflguratioo 

Laminates loaded to fracture in 

Number of ply 
stress states 
at fracture 

TensloD 

Siear 

Compression 

X 

y 

XY 

X 

Y 

Per 

laminate 

Total per 
system 

^Graphite/epoay ASV3501 (A&^E) 








(Ogl 

x 

X 

X 

X 

X 

5 


[0/±4V(V±45lg 

X 

X 

X 

X 

X 

15 

60 

I Oy'i 45/02,^90/01 g 

X 

X 

X 

X 

X 

20 


|(O/*4V90>2|g 

X 

X 

X 

X 

X 

20 


« 

^Boroo/epo]^' 5505/4(B/E) 









X 

X 

X 

X 

X 

5 

25 

(G/*4S/90)g 

X 

X 

X 

X 

X 

20 


^E-Glaas/epoxy <B-G/E) 








{Ogl 

X 

X 

X 

X 

X 

12 


1 ^ 30)31 5 

X 

X 

* 

X 

X 

20 


[ (* 45>2 1 jj 

X 

X 

X 

X ^ 

X 

40 

142 

M*60)2)g 

X 

X 1 

X 

X 

X 

30 


[SOgl 

X 

X 1 

X 

X 

X 

40 


^Graphlte/epoxy MOO-l/ERLA 








4617 MOXLI/E 









X 

X 

X 

X 

X 

5 


{O/*60j g 

X 





21 


( 0/±43/90) g 

X 
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^Flat specimens subjected to uniaxial load. 

"Thin tubes subjected to uniaxial and comUned loatk mith icplicates). 

specimens subjected to uniaxial loads at 0°, IS'^. 10®, 45®, 60®, 75® 
and 90® to 0-ply. 


TABLE IIL - COMPARATIVE VALUES OF THE Aj COEFFICIENTS IN EQ. (1) FOR THREE 


DIFFERENT EVALUATION CONDITIONS 


Evaluation 

Coefficients normalized with respect to A^ 

condition 

Af/Ai 

A^/Aj 

A 3 /A 1 

A/A, 

A 5 /A, 

Ag/Ai 

A 7 /A, 

kPa/(ksi) 

Ag/Ai 

kPa^ (ksl) 

Vh 

kPa/ (ksi) 

Aj>< 0 

1.0 

51.8 

407 

- 10.0 

1111 

-926 

-281 

2551 

460 







i-‘i0,7) 

(370) 

( 66 . 7) 

A5=\ = ® 

1.0 

51.8 

407 

- 10.0 

0 

0 

-281 

2551 

434 







i-40.7) 

(370) 

(62.9) 

A 5 “ Ag == A^ = 0 

1.0 

51.8 

407 

- 10.0 

0 

0 

-281 

2551 

0 







(-40.7) 

(370) 




*Aj = 0.57x10"® MPa'* (0.27x10"'* (ksl)"^). 
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Figure 1. Flow chart for a prxedure to deterrtjne the n situ ply strengtt.s 
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